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Abstract 

We study the coupled Einstein-Klein-Gordon equations for a 
complex scalar field with and without a quartic self-interaction 
in a curvatureless Friedman-Lemaitre Universe. The equa- 
tions can be written as a set of four coupled first order non- 
linear differential equations, for which we establish the phase 
portrait for the time evolution of the scalar field. To that pur- 
pose we find the singular points of the differential equations 
lying in the finite region and at infinity of the phase space 
and study the corresponding asymptotic behavior of the so- 
lutions. This knowledge is of relevance, since it provides the 
initial conditions which are needed to solve numerically the 
differential equations. For some singular points lying at in- 
finity we recover the expected emergence of an inflationary 
stage. 
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1. INTRODUCTION 



The recent developments in particle physics and cosmology suggest 
that scalar fields may have played an important role in the evolution 
of the early Universe, for instance in primordial phase transitions, and 
that they may constitute part of the dark matter. Moreover, scalar 
fields are predicted by most of the particle physics models based on 
the unification of the fundamental forces, as for instance in superstring 
theories. Scalar particles are needed in cosmological models based on 
infiation, whose relevance is supported by the results of the COBE- 
DMR measurements that are consistent with an Harrison-Zel'dovich 
(scale invariant n=l) spectrum |1|]. These facts, in particular infia- 
tion, motivated the study of the coupled Einstein-scalar field equations 
to determine the time evolution and also the gravitational equilibrium 
configurations of the scalar fields. The latter one in particular for mas- 
sive complex scalar fields, which may form so-called boson stars 0, ^ . 

A detailed study of the solutions of the Einstein equations for a 
homogeneous isotropic Friedman-Lemaitre Universe with a real scalar 
field has been done in particular by Belinsky et al. 0, ^, H, see also 
ref. [[7[. In this paper we extend, following ref. |^, ^, these investiga- 
tions to a complex scalar field. 

This analysis is important in order to see the degree of generality 
of solutions possessing an infiationary stage and also due to the fact 
that these solutions constitute the background, starting from which 
one can study in the early Universe the time evolution of perturbations 
for the scalar field H, A fact this, which is of relevance if scalar 
fields make up part of the dark matter. If this is the case they may 
also form compact objects, such as bose stars, or trigger the formation 
of observed large scale structures in the Universe. 

Since we consider a complex scalar field with or without a quartic 
self-interaction in a curvatureless Friedman-Lemaitre Universe, we get 
for the Einstein-Klein Gordon equations a set of four first order non- 
linear differential equations, for which we study the phase portrait for 
the time evolution of the scalar field. We first determine the singular 
points of the differential equations and then find analytically the asymp- 
totic behavior for the solutions nearby these points. For the singular 
point lying in a finite region of the phase space, we can use for m 7^ 0, 
in an adapted coordinate system, the averaging method in order to get 
the asymptotic behavior of the solution, whereas, for the points lying 
at infinity, we first compactify the phase space on the lower hemisphere 
of the 3-dimensional sphere. We can then apply the Poincare-Dulac 
theorem to get the corresponding asymptotic behavior. From the so- 
lution we then see if there is infiation and how long it lasts. For some 
singular points lying at infinity, we recover the expected emergence of 
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an inflationary stage. 

The paper is organized as follows: in section 2 we present the basic 
equations which we will use. In section 3 we first study for a massive 
scalar field the singular point lying in the finite region of the phase space 
and then the singular points lying at infinity both with and without 
a quartic self-interaction. Section 4 is devoted to the massless scalar 
field with a quartic self-interaction and a short summary concludes the 
paper. 



2. BASIC EQUATIONS 



We consider a massive complex scalar field with quartic self-interaction 
in a Friedman-Lemaitre Universe with the following action 

^ = / (-TTT?^ + ^z^^^''^* + + A(w*)')v^ d'x , (1) 

J iovrG 

where g is the determinant of the metric 

ds^ = g^jf'e" = -e°e" + s^jO'e^ , (2) 

with 



^ ^ 4' J i=i 



and is the dual basis of 9^''. By varying the action with respect to 
g^^ we get the Einstein field equation 

G^, = SttGT^, , (3) 

with 

T^u = e^v5ei,v5* + e,.v5e;,v5*-5'^,.(5'"'^eaV5e/3V5* + m^V9V9* + A(v9V3*)^) . (4) 
The (00) component of eq.(H) leads to the constraint equation 



H' + -, = ^(W* + ^'VV* + Kw*f) , (5) 
a o 



with H = d/a and dot means derivative with respect to time. For the 
(ij) component we have 

{-2H - 3H^ - ^)5ij = SnGiifiP* - m^ipcp* - \{ipip*)^)5^j . (6) 

By varying the action with respect to and we get the Klein-Gordon 
equation 

(f + 3Hip + m'^ip + 2X{ipip*)ip = (}, (7) 
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and its complex conjugate. The system is fully determined by the 
independent eqs.(|]) and (^. In fact, one can easily show that eq.(|) 
follows from eqs.(^ and (|^. For the massive scalar field case we use 
the following dimensionless variables 



t ^ 7] = mt , 

\ V A — ^ 

^ ^ ■'^ — 87rGm2 ' 



_ m ' 
h h — ^ 



This way, we get for eqs.(|D and the following set 
k 



a 



2+^2 = yf + yl + xl + xl + 3A(x? + xlf , (9) 

yi = — 3zyi — xi — 6Axi{xl + xl) , (10) 

x'l = yi , (11) 

y'2 = - 3zy2 - X2 - QAx2{xl + xl) , (12) 

x'2 = y2 , (13) 



where prime means derivative with respect to rj. The only singular 
point (defined as the point {xi, X2,yi, 1/2) for which the right hand side 
of eqs.(P) - (|13D vanishes), which we denote by A, lying in a finite region 
of the phase space (defined by ip, (p or equivalently Xi, X2, yi, 1/2) is the 
coordinate origin. Eqs.(P) - ([l3[) are invariant under the transforma- 
tions 

a) Xi ^ — Xi and yi ^ — yi , 

b) X2^ - X2 and y2 ^ - 1/2 , (14) 

c) X2 ^ xi and 1/2 ^ yi ■ 



For every solution which describes an expanding Universe (i.e. z > 0) 
there is a corresponding solution describing a collapsing Universe. This 
can be seen by performing one of the following transformations on the 
set of eqs.(ig) - (|T§) 



T] - 


y -7] 


z 


► —z 


Xi - 


-Xi 


X2 - 


-X2 



V - 


> -T] 






> -7] 


z — 


y —z 






► —z 












Xi - 


-Xi 




1 


-yi 


y2 - 


-1/2 




[ X2 - 


-X2 



b) { ^ . \ , c) ^. , rf) ^. .(15) 
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For the massless scalar field case we use the dimensionless variables 



^0 



SttG 



V2 ^ a;io + ix2o = y^V^ , 

^ yio + iy2o = ^0 , (16) 
H -^zo = VS^H 
k ko = SttG k , 



and A remains unchanged. This way, we get for eqs.(^ and (|^) the set 

^ + ^0 = Z/io + 1/20 + 3A(x?o + a;^o)% (17) 

y'lo = - 'izoViQ - Q\xw{xIq + X20) , (18) 

a^'io = yw , (19) 

y'20 = - 32:02/20 - QXx2oixjQ + xIq) , (20) 

3^20 = , (21) 

where prime means here derivative with respect to rjQ. These equations 
are also invariant under the transformations given by eqs. Ql^) and (pTSl), 
and the coordinate origin is the only singular point lying in the finite 
region of the phase space. 



3. MASSIVE SCALAR FIELD IN A CURVATURELESS FRIEDMAN-LEMAI- 
TRE UNIVERSE 

As next we study the asymptotic behavior of the solutions of eqs. (y)- 
(p^) nearby the singular points. Due to the increasing complexity of 
the analysis involved for a complex scalar field with respect to a real 
one, we restrict ourselves to the case k = 0. In section 5 we briefly 
comment on the extension to = ±1. With = in eq.(H) we can 
then rewrite eqs.(|T0|)-(|13]) in spherical coordinates, defined as follows 



xi = r cos-i^a cos^^i , 
yi = rcos'ds sin-i^i , 
X2 = r sin cos '§2 , 
1/2 = r sin sin '62 , 



[22) 



with i?!, 192 G [0, 2t[) and t^s G [0, vr). This way we obtain 

= + 3Ar^(cos2'(93Cos2^9i + sin2'(93Cos^^?2)% (23) 
i)i = —3z sin t^i costal — 1 — 6Ar^(cos^ -i^a cos^ "i?! + sin^ "i^a cos^ 192) cos^ -(9(2,4) 
= — 32; sin'(92 cos'(92 — 1 — 6Ar^(cos^ "i^s cos^ -i?! + sin^ t^s cos^ -(92) cos^ ?9^2,5) 
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-d^ = sin'i93Cos'i93[— 32;(sm^'i92 — sin^-i?!) 

— 6 Ar ^ (cos^ cos^ + sin^ cos^ ^2 ) (sin '(92 cos ^2 — sin cos 'i9{!^6) 
r = — 3rz(cos^ -i^s sin^ "i^i + sin^ -i^s sin^ 'i?2) 

— 6Ar^(cos^ -i^s cos^ -i^i + sin^ ^3 cos^ -(92) x 

(cos 'di sin -i^i cos^ + cos {^2 sin "(92 sin^ -i^s) . (27) 



Due to the transformations given in eq.(p!4D which leave the equations 
invariant, we can restrict the domain of definition for t^i, ■}}2 and t^s 
respectively to [0,7r), [0, vr) and f,vr^. At the singular point A lying 

at r = the above equations reduce to = —1, ■t}^ = —I, "^'^ = and 
r =0. To get the asymptotic behavior of the solution nearby A we 
apply the averaging method (for details see ref. fT^)- Thus, we have 
to solve the differential equation 

r' = -lr\ (28) 

which is obtained by averaging eq.(P7D over the angular variables. We 
get the asymptotic behavior for the solution nearby A 



xi = cosT^so cos(r7 - 771) , 

yi= - cost^so sin(77 - ?7i) 

X2 = ^ sin t^so cos(r7 - 772) , 

2/2 = - 1^ sin ??3o sin(r7 - 172) 



(29) 



where ■(^so, rji and 772 are integration constants. We see that A is an 
asymptotically stable winding point. Setting 7^30 = corresponds to 
consider only a real scalar field and we recover the solution discussed 
by Belinsky et al. in ref. 0, |^, |^. 

We study now all singular points lying at infinity in phase space. 
First, we consider the case with no quartic self- interaction term. 

3.1. Properties of the singular points lying at infinity for A = 

In order to find the singular points lying at infinity we perform a 
transformation which maps them on the boundary of a unit 3-sphere, 
defined as follows 



dr] 



with p G [0,1). With this transformation eqs.(|5) - (127|) become 



(30) 



^ = -Sp^il- p){cos'^^3sm'^^i + sm^'&3sm'^^2) , (31) 
dr 
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dr 

dA 
dr 
dd 



dr 



-- -3psin?9icos^9i - (1 -p) , (32) 
= — 3psin't92 cos'i92 — (1 — p) , (33) 
^ — 3psin'i93 cos793(sin^'(9i — sin^'i92) ■ (34) 



Since the system of difFerential equations is well defined for p = 1, we 
can extend the domain of definition for p to the boundary. This way, 
we have now a compactified phase space. At infinity (i.e. setting p = 1 
in the above eqs.(|3l|)-(|3^) we find two singular curves and two singular 
points, which are given by 

, ^ < < vr , 
7T vr 

To study the asymptotic behavior of the solutions nearby li is rather 
involved, due to the fact that every point lo = (1, 0, 0, t^so) G h is non 



h : 




= 0,^92 








k : 


^1 






^1 




Pi ■■ 




P2 ■■ 


^1 


= 0,^2 



hyperbolic (see ref. |Tl| for details). We first perform a variable shift 
in eqs.([3l|)- (|3^) defined as follows: 5p = p — 1, 6'^i = -di, 8'd2 = '^2 
and 6-})^ = "^3 — "^so- We then expand the differential equations in the 
new coordinates up to third order in a sufficiently small neighbourhood 
around the singular point Iq. Next, we define the linear coordinate 
transformation: x = 6p, y = —^Sp+S-di, v = —^5p+5'd2 and w = S-ds, 
such that we get a set of differential equations of the form 

{y)'=Dy + p, (35) 

where y = (x, y, v,w). D is a diagonal matrix and p is a vector whose 
components are polynomes in x, y, v and w containing monomes of 
degree 2 and 3. We can now use the Poincare-Dulac theorem [12| to 



classify the monomes in p into resonant and non-resonant ones. For 
the non-resonant monomes of degree n, there is a polynomial change of 
coordinates of degree n, so that they are transformed into polynomes 
of at least degree n + 1. This is not the case for the resonant monomes. 
The polynomial change of coordinates is found by solving the so-called 



homological equation (see ref. |jT2| for more details). Performing the 
polynomial change of coordinates on the differential equations, the non- 
resonant terms of degree 2 and 3 become of higher order and are thus 
neglected. Retaining only terms up to third order we obtain the set of 
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equations 



3^ = - 3y - 3£y + ^sin^?93o(x^y - , (37) 
ar 3 

^ = — d,v — 3x5 cos^ 'd^Q^x'^y — x^v) , (38) 

(ir 3 

dw 

dr 



, (39) 



where {x,y,v,w) = {x,y,v,w)+ polynomes of degree 2 or higher in 
{x,y,v,w). Since the solution of eg. (|36|) does monotonically increase 
as a function of r, Iq is a saddle point. The above equations can now be 
solved analytically and furthermore by transforming back to the original 
variables one gets for the asymptotic behavior nearby any point Iq of h 

f = ^ . (40) 
^ = ^ e-" , (41) 

H = ^. (42) 

for t — »• — oo, corresponding to the initial singularity and where 
Ma = /} ^ - This solution corresponds to an outgoing line, the so- 

^ vottG- 

called separatrix. Using eq. (^) and the definition of H we have 



exp -:m'(tf - ti) (43) 



for tf > ti. 

The effective equation of state near Iq tends to e = —p, where e = Tqq 
and p = I (Til + T22 + T33). The solutions near the separatrix are 
characterized by the fact that 00* -C m^ipip* and the phase ^^so being 
constant. If for a time tf—ti a trajectory T satisfies these two conditions 
and, at time t/, T is near to the separatrix, then it follows with eq.(^) 

that if ~ ^^^^ He^"^^". Using eq.(^ it is easy to show that 



a{ti) 



<f{ti) + Ci 



Cfitf) + Cl 



(44) 



for tf > ti, with Ci = ^^^^ e"^^°. Every trajectory which lies close 
enough to the separatrix will thus meet the criteria for inflation. 
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To establish the asymptotic behavior of the solution nearby the sin- 
gular saddle point pi, the same strategy as for the line li has to be 
applied. We therefore give here only the result, which is 

V = (45) 
^ = ^ , (46) 

* = ^ • '^^) 

for t ^ —oo. This solution corresponds also to an outgoing separatrix. 
A more detailed analysis shows that the singular point can only be 
reached if "i^s = f . Otherwise, starting from points nearby pi with 
"i^s 7^ 7r/2 the phase of the scalar field varies strongly. In this region of 
the phase space inflation is driven by the imaginary part of the scalar 
field. Also here the phase of ip remains constant along the separatrix 
and we get the same equation of state as for the previous case. Applying 
the transformations defined in eq. ([T^ ) on eq. (|45|) , we get instead of pi 
the singular point pi = (1, 0, |, 0), for which the corresponding solution 
is now real. Hence, the analysis made in ref.0, |^, ^ applies here as well. 

For all points b = (1, |, |, i}^) in I2 we can directly solve the linearized 
differential equations. It turns out that remains constant and that 
on the plane i)^ = "i^so = const the solution expands away from b. 
Using the inverse transformations of eqs.(pUD and (^, we obtain the 
asymptotic behavior of the scalar field and of the Hubble parameter 

.^^.„(i)e-. (48) 

for t ^ 0^ corresponding to the initial cosmological singularity and 
where to is an integration constant. The equation of state nearby these 
points tends to e = p (i.e. stiff matter). From eqs.(||) and (H) we get, 
as long as ipip* ^ m'^ipip*, that H = —3H^. Solving this differential 
equations and using the definition of H, we obtain as expected 

^M=(^)*. (51) 

a{ti) \tiJ 

Following the same method used for the points in I2, we get for the 
saddle point p2 



(52) 
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H 



-C2mH i 



1 

3t 



(53) 
(54) 



for t — s> 0+, with C2 and to being integration constants. The results of 
the analysis for I2 do also apply here. 

This completes the study of the phase portraits for A = 0, for which 
we found two singular curves /i, I2 and two singular points pi, p2- AH 
other singular curves and points, due to the transformations given in 



eq.(|T4D, can be reduced to one of this cases. For the curve li and the 
point pi the solutions correspond to an outgoing separatrix where infla- 
tion occurs. Along these separatrices the phase of if remains constant. 
For the curve /2, for which the effective equation of state corresponds 
to stiff matter, the phase remains also constant. Setting = in 
the solutions around li and 12-, we recover the results found by Belinsky 
et al.[^, ^ ^ for a real scalar field. Contrary to the previous ones the 
solution around p2, for which we also get an effective equation of state 
for stiff matter, can not be obtained by simply adding a constant phase 
to the corresponding asymptotic solution for the real scalar field. As 
next we turn to the case where there is a quartic self-interaction term. 

3.2. Properties of the singular points lying at infinity for A 7^ 

We perform on eqs.(p^)-(|27|) the transformation defined by eq. (pO|) . 
We then obtain the set of equations 

= — 3p^/(cos^ i^s sin^ -(^i + sin^ t^a sin^ ??2) 

dr 

„2„Q I "'-^^2 cos'i92 sin^ t^s) , (55) 

,(56) 



d^i 




1 


dr 


1 


-P 


d^2 




1 


dr 


1 


-P 


d^3 




1 


dr 


1 


-P 



— QKp^g (sin -(^i cos t?! cos^ 'd^, + sin '82 

{1 — pY + 3pf sin cos -(^i + GAp^g cos^ 

{1 - p)'^ + 3pf sin ^2 cos ^2 + QAp'^g cos^ ^2] ,(57) 
3pf sin cos -i^a (sin^ '§2 — sin^ -§1) 



+ 



where 



and 



/ = V'(l-p)2 + 3ApV 
g = cos^ 'ds cos^ 'di + sin^ t^s cos^ '^2 ■ 



(59) 
(60) 
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Since the right-hand side of eq. (|55| ) is well defined and continuous for 
p = 1, we can find all sets of angles ("i^io, ^^20; ''^so) for which the condition 
^ = at p = 1 is fulfilled. For each solution {-^lo, 1^20, ^30) we have to 
check if linip^i- fi{p, ^10, ^20, ^30) = 0, where fi stands for the right- 
hand side of eqs.(|56D - (|58D. If this is the case, the point (1, -dio, {^20^ ^^so) 
is a singular point at infinity. Again, we find two singular curves and 
two singular points 



Li : 


^1 


= arctan 




P - 




L2 : 


^1 




Pi : 


^1 


TT 


P2 : 


^1 


= 0,!)2 = 



-2V3A\ „ /-2V3A\ TT 

+ 71 ,V2 = arctan ( | + tt , — < vs 



2 



TT TT 

-,-<»3<7r,p^l, 

/-2V3A\ , TT 

arctan I j+vr, ^$'3 = -,p-^l, 

^ Q ^ 1 

-,^3 = ^,P^l. 

Notice that the limit A ^ is rather subtle, since in the differential 
equations there are terms involving A/(l — p), which are not defined 
when p — i> 1 and A — > simultaneously. We expand the differential 
equations around an arbitrary point Pq of Li defined by the following 
coordinates 



/-2V3A\ , /-2V3A\ 
p, Vi = arctan + vr, -1/2 = arctan + vr, = V30 , 



where '^^q is in [|, vr). In order to have finite partial derivatives, when 
p tends to 1, the angular variables 'di and '^2 are kept fixed. When 
linearizing eq.(^8|) around Pq, we see that also the angular variable "i^s 
has to remain constant. Inserting the values of the angular variables 
defined in eq. (p2D , we get for the scalar field 

-3Mpm e^'^so 
if = I = r , (65) 

'4AM2 + 3m2 




ip = -2^-M,ip, (66) 
for r —>■ 00, where r is a dimensionless parameter defined by eq. ([2^) . 



Since the phase of (p is constant, we can integrate eq. (|66D and thus get 
as a function of t rather than r. This way we get 

<^ = (^oe'''^«exp |-2Mp7|t| , (67) 
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where t — > — oo and ipo is a negative integration constant. Using eq.(|^), 
one gets the asymptotic behavior for the Hubble parameter. This so- 
lution corresponds to an outgoing separatrix and the equation of state 
tends to e = —p. A trajectory, which lies sufficiently close to the sep- 



aratrix, can be parametrized by ~ a(p, where a ~ — 2Mp-y/A/3. As 
long as ~ X^ipip*)"^ , with eq. (|^ we find along the trajectory 

'■"^>-expf „ -f^, (t,-tA. (68) 



for tf > ti. The factor in the exponential multiplying the time difference 
is just the Hubble expansion rate, which is positive and tends to infinity 



as a reduces to — 2MpyA/3. Thus, any solution which lies sufficiently 
close to the separatrix will go through an infiationary stage. 
Using the same strategy as before for the point Pi, we obtain 



^ = i^oexp{-2MpJ-t] , (69) 




where t —oo and ipo is a negative integration constant. This solution 
corresponds also to an outgoing separatrix. The singular point can only 
be reached if "i^s = |, thus the real part of the scalar field vanishes. The 
analysis made for Li applies as well, again there is an infiationary stage 
if the solution gets close to the separatrix. 

For the lines L2 we expand the set of differential equations to first 
order around a point B defined by the coordinates (p, |, |, -i^so). For 
p —>■ 1 the linearized equations reduce to those obtained for the line l2- 
Therefore, the solutions are given by eqs.(^) - (|50D. The inclusion of 
a quartic term in the lagrangean does not effect the singular line. This 
is expected, since near L2 the potential term is negligible compared to 
the kinetic term 00*. The same remark holds for the singular point P2, 
so that the asymptotic behavior is given by eqs.(|52|)-(p^. 

We see that the presence of a quartic self-interaction term does not 
change substantially the main features of the phase portrait. We find 
again two singular curves and two singular points. All other solutions 
can be reduced to these by the transformations given in eq. (p!4D . The 
solutions which correspond to outgoing separatrices show an infiation- 
ary stage. The behavior of the solutions around these separatrices is 
not much affected by the quartic term. Its main influence is to shift 
the position of the singular points, where inflation occurs. For the solu- 
tions for which no inflation occurs, that is around L2 and P2, the results 
are the same as for li and p2 respectively. We obtain the solutions of 
Belinsky el al. in ref. for the real scalar fleld by setting ■(^a = as 
mentioned in section 3.1. 
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In Fig.l and 2 we plot the numerical solutions of the differential 
equations (^^) - (|13|) for different values of A. For all solutions we take 
the same initial conditions as given by the asymptotic behavior near 
the singular point p2- The plots are valid for every value of m 7^ 0, since 
eqs.(0) - (|T3|) do not depend on it. The numerical solutions have only 
a physical meaning in the region where Tqq < M^. Hence, depending on 
the value of m, it might be that only a part of the Figure is of physical 
relevance. 

As next, we analyse the asymptotic behavior of the solutions of 
eqs.(p!^)-(|2T|) for a massless scalar field with a quartic self-interaction 
nearby the singular points. 



4. MASSLESS SCALAR FIELD IN A CURVATURELESS FRIEDMAN-LE- 
MAITRE UNIVERSE 

Again the coordinate origin is the only singular point lying in the 
finite region of the phase space. This point is asymptotically stable. 



since we have the following Liapunov function (see ref. |13| for details) 

Kxio, X20, yio, 2/20) = yfo + Vw + 3A(2;^o + xj^y , 

for which the derivative with respect to r/o is strictly negative, except at 
the origin. Applying the transformation defined in eq.(p2D on eqs.(|l^)- 
(0), we get 



12 



X 



^1 



|l + 3Ar^(cos^ -i^s cos^ -^i + sin^ i!}^ cos^ ^^2) - 
(cos^ ^3 cos^ ?9i + sin^ cos^ ^^2) } , (70) 
— 82:0 sin-i^i cos-i^i — 6Ar^(cos^ 'ds cos^ -i^i + sin^ -ds cos^ 'd2) cos^ 'di — sin^ ^1) 
—3zo sm'd2 cos'i92 — 6Ar^(cos^ cos^ i^i + sin^ i^a cos^ ■i?2) cos^ '(92 — sin^ ^2,) 
sin-i^s cos-i^sf— 32;o(sm^ ^^2 — sin^^^i) 



— 6Ar^(cos^ t^s cos^ i?! + sin^ cos^ 192) (sin ^92 cos'(92 — sint^i cost^i) 
+ (sin {^2 cos '(92 — sin -i^i cos i^i)] , 

r = — 3rzo(cos^ "i^s sin^ i?! + sin^ t^s sin^ -(92) 

+ r(cos'(9i sin-i^i cos^^^a + cos'(92 sin'i92 sin^^^a) 

— 6Ar^(cos^ -i^s cos^ + sin^ 7^3 cos^ 7^2) x 
(cos sin 'di cos^ -(^a + cos '^2 sin -(92 sin^ -(^s) . 

In a sufficiently small neighbourhood around the origin, which we de- 
note by W, we consider the projection of the solutions of eqs.([7lD-([74D 
on the (xio, |/io)-plane. The angular variable on this plane is and 
its behavior is given by the solution of the equation = — sin^t^i, 
which is just eq. ([TTD with r = 0. Therefore, knowing that the coordi- 
nate origin is asymptotically stable and that i)i is almost everywhere 
strictly negative, it follows that the solutions are winding towards the 
point Aq = {xio = 0, yio = 0). One gets the same behavior when the 
solutions of eqs.([7Tl)-([7^) in W are projected on the (x2o, ?/2o)-plane. 

We now turn to the singular points lying at infinity. One has to 
apply on eqs.(^)-(74) the transformations given in eq.(0), where f] 
has to be replaced by tiq. We obtain the following equations 



dp 
dr 



JL = -3p^f{cos^^3sm^^^ + sm^^ssm^^2) 



+ 



-6Xp^g + p{l-p)' 



sin 'di cos cos^ + sin i!}2 cos 'd2 sin^ 



(75) 



(Ml 

d^ 
dr 

d^ 
dr 



1 



1-p 
1 



1-p 

cos -i^s sin -i^a 



1 



3pf sin cos -i^i + QXp^ g cos^ + (1 — p)^ sin^ ^{''^Q) 
3pf sin 7^2 cos^92 + 6Ap2 g cos^ ^^2 + (1 - P? sin^ 
[3pf{sm'^^2-sm\h) 



P 



+ 



6Ap^ — (1 — p)^ (cos -(92 sin -(92 — cos?9i sin-i^i)! 



(78) 



(73) 



(74) 
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where from now on 



/= v/(l-p)^(l-^7) + 3A/^V (79) 



and g is still given by eg. (|60|) . The singular points lying at infinity of 
the phase space are found using the same method as for A 7^ 0. We get 
again two singular curves and two singular points, which we denote by 
^10; ^205 Pio and P20 and their coordinates are given by eqs.(pTD-(|5^), but 
now with A replaced by A. To obtain the asymptotic behavior around 
these singular points we apply the same method used in section 3.2. 

The asymptotic behavior around a singular point lying in lio is given 
by eqs.(p3D-(pBD, but where m is now replaced by Mp. The angular 
variables must be kept fixed in order to have finite partial derivatives, 
when p —>■ 1. The analysis made for Li is also valid here and /lo 
has infiationary stages. Setting -^3 = we obtain automatically the 
asymptotic behavior for a massless real scalar field with a quartic self- 
interaction for which we recover the infiationary stage. A fact this 
which was established heuristically by Linde in ref. . 



One gets the behavior of the solutions around the singular point piQ 
from the one near Pi in the same way as discussed above for /lo from 

The asymptotic behavior of the solutions around all points c = 
(1, |, |,'i93o) of I20 or around the point P20 is found directly by solv- 
ing the corresponding linearized differential equations. For the singular 
points on /20 it turns out that the solutions are given by eqs.(^)-(^), 
whereas for p2o we get the solutions from eqs. (p2|) -(^4|) inserting m = 0. 



The form of the solutions around the line /20 and P20 was expected to 
be similar to that obtained for I2 and p2, because the potential term is 
negligible with respect to the kinetic energy term ip<f*. 



5. CONCLUDING REMARKS 



The extension of the above analysis to the singular points for = ±1, 
although in principle straightforward, is much more involved. One 
could, for instance, using eq.(^ eliminate the curvature term in eq.(|) 
and consider this modified equation together with eq.(0). This gives 
then a set of five non-linear first order differential equations. As a con- 
sequence, when performing the transformation to spherical coordinates 
needed in order to compactify the phase space, one gets one additional 
angular variable. It turns out that around some singular points the ex- 
pansion of the differential equations must be done at least up to fourth 
order. The only singular point lying in the finite region of the phase 
space is at the coordinate origin. We conjecture that the asymptotic 
behavior near the singular points lying at infinity with k ^ will not 
fulfill the criteria for infiation. This does not imply that infiation can 
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not occur in an open or a closed universe, but that every trajectory 
must come close enough to one of the separatrix found for = in 
order to go through an inflationary stage. This fact has been shown for 
the real scalar field case (see ref. [§|, |, H). 

In this paper we have extended to complex scalar fields the analysis 
of the initial conditions in an homogeneous and isotropic Friedman- 
Lemaitre Universe. The main features found for real scalar fields hold 
also for complex scalar fields, in particular the existence of infiationary 
stages. The fact that along the separatrices the phase of (p remains 
constant is important and shows that infiation is essentially driven by 
one component of the complex scalar field. Therefore, the results on 
infiation valid for a real scalar field (see for instance ref.|T5| and refer- 



ences therein) apply also on the component of the complex fields which 
drives infiation. The behavior around the singular points p2i P20 and 
P2 is more involved and can not be obtained by just adding a phase to 
the corresponding solutions for the real scalar field. We also notice that 
for a massive scalar field the presence of a quartic self-interaction term 
does not change substantially the main features of the phase portrait. 
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Figure Captions 

Fig. 1: Plot of the dimensionless variable yi as function of Xi for 
different values of A. It corresponds to the phase portrait 
of the real part of the scalar field. All curves start with the 
same initial condition given by the asymptotic behavior near 
P2 (see eq. (^2])). The real part of the scalar fields depends 
on the value of C2, for which we take the value C2 = 1. 

Fig. 2: Plot of the phase portrait of the imaginary part of the scalar 
field in dimensionless variables for the same cases as in Fig.l 
The asymptotic behavior of the imaginary part near p2 de- 
pends on the value of to, for which we choose to = 100. 
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